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Abstract 

Starting from the linear flow of homogeneous fluid, five modes are defined as eigenvectors of the 
basic system of conservation laws. Quasi-plane geometry is considered. Projectors that separate overall 
perturbation of the fluid into specific modes are calculated and applied to nonlinear flow. Dynamic 
equations for the interacting modes are obtained in the frames of the method. A particular case 
of streaming caused by acoustic pulse is considered, illustrations on temporal behavior of streaming 
velocity and streamlines are presented. 

PACS numbers: 43.25.Nm 
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INTRODUCTION 



Streaming is known as a bulk vortical movement of the fluid following an 
intense acoustic wave. Streaming appears in closed and open vessels in two- 
or three-dimensional geometry and is observed in viscous fluids only. Though 
vortical movement may exist in the linear flow, the reason of streaming are 
the nonlinear losses of momentum of acoustic field. A reader is referred to the 
comprehensive reviews on this subject (if, J2|, J3|. Most recent investigations 
)oth theoretical and experimental relate to steady (quasi-steady) streaming 
It is well-known that streaming is a slow process in comparison to the 
ultrasound and may be separated from the originating sound by time-averaging 
over integer number ofperiods of sound. The averaging is a starting point 
of the modern theory |2|, j^]. A driving force of streaming is quadratic and 
therefore does not vanish during the temporal averaging. Some inconsistencies 
of the approach based on temporal averaging are obvious. It was first pointed 
out by Rudenko and Soluyan that studying acoustic streaming while assuming 
incompressibility was inadequate J2|. The effect of compressibility of the fluid 



has not been discussed in depth. Recently, the importance of compressibility 



has been demonstrated for flow over plane rigid boundary [5] . It is more evident 
in gases and leads to larger streaming velocities. 

We may add that incompressibility neglects not only acoustic waves but also 
another slow compound of the overall flow, namely the entropy modal field that 
appears in one-dimensional geometry of flow even. Entropy mode appears due 
to losses of acoustic energy and therefore has a different nature than streaming. 
Velocity of this field is small in comparison to streaming JSj, the main feature 
of entropy field is slow isobaric growth of temperature that leads to the new 



background of acoustic wave propagation. This effect called by acoustic wave 
self-action is well-studied (l6( . Excluding the entropy mode needs a neglecting 
of equation of energy balance that reduces an initial system of conservation laws 
and allows to consider flows with unperturbed density. At last, an excluding of 
the energy balance does not allow to account heat conductivity though it has 
been proved that the effect of the heat conductivity could not be discarded in 
a study of temperature variation associated with the streaming p. Actually, 
this approximation is well understood in a typical liquid like water but should 
be revised for other liquids |7|. 

The dynamic equations are derived by proper combining the conservation 
laws in differential form. These calculations may be considered as a certain 
type of projecting at the chosen type of motion. We propose to determine 
modes and correspondent matrix projectors from the very beginning at the 
level of linearized system of conservation laws. It seems reasonable to start 
from the full system of conservation laws including the energy balance equation 
and the continuity equation to avoid possible inconsistencies. 

Though ultrasound of many periods is most efficient source of streaming, any 
other acoustic source in thermoviscous flow gives rise to other modal fields. The 
modern theory fails to treat non-periodic sound. That is one more reason to 
find new methods (analytical more desired) in fluid dynamics and particulary 
in the theory of streaming. The idea to decompose the linear flow into specific 
modes is not novel and has been exploited for a long time, see the paper by Chu, 
and Kovasznay [9] and referred there papers, where homogeneous background 
with sources of heat and mass are considered. This paper is almost the only 
one in which interaction between modes is introduced locally (see also report by 
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Kaner, Rudenko, Khokhlov [10], where a first step has been made in a manner 
to introduce interaction of opposite acoustic waves via connection equations). 
The account of nonlinear interactions by Chu, and Kovasznay (2| was performed 
by the simple perturbation theory - by decomposition of each variable in series 
of small parameters of nonlinearity and dissipation. 

The concrete ideas of doing it automatically by using of projectors in the 
wide variety of problems was realized for flow over inhomogeneous media like 
bubbly liquid Jll| , flows affected by external forces including the gravitational 
one which changes the background density and pressure |l2j |. both one- and 
multi-dimensional problems [13]). The principal advance is an expansion of the 
ideas into area of nonlinear flow: to get nonlinear coupled evolution equations 
for interacting modes and to solve the system approximately. The choice of 
the subspaces by the projecting (fixed links between variables) is the principal 
point of evolution equations derivation. It results in the expansion of projection 
of the evolution operator as a function of a small parameter , that is equivalent 
to the nonsingular perturbation theory l^j]. Any weak interactions of modes 
in fluid dynamics may be considered. Streaming is a particular case of the 
interaction of modes when the acoustic mode is dominant and inverse influence 
of the growing vortical end entropy modes neglected. 
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MODES AND PROJECTORS OF THE QUASI-PLANE THERMOVISCOUS FLOW. 
DYNAMIC EQUATIONS OF LINEAR AND NONLINEAR FLOWS. 



The mass, momentum and energy conservation laws for the thermoviscous 
flow are the starting point: 



dt 

dv 



§ + v(H = o 



§ + (uV)i?J = -Vp + r]Av+ (q + f ) V(W) 
+ (vV)e] +pW- X AT = q (W) 2 + I (H + || - f^gg) 2 



P 

p 

In this system p,v : p,e, T denote pressure, velocity, density, internal energy 
per unit mass and temperature relatively, q, 7/, x are bulk, shear viscosities and 
thermal conductivity (all supposed to be constants), X{ - space coordinates, t- 
time. A system ( TTI) should be completed by two thermodynamic relations for 
equilibrium processes in fluid, e{p, p),T(p, p) . For simplicity, an ideal gas will 
be considered treated by the following relations: 

p p 

e = ~^Ty T = p( 7 - Wv ( ' 2) 

with 7 = Cp/Cv' } Cy,Cp being heat capacities per unit mass under constant 
volume and pressure correspondingly. Any other fluid may be considered as well 
by expansion of thermodynamic relations into the Taylor series in the vicinity 



of equilibrium state (22] 



The quasi-plane flow along y-axis will be considered. That allows to introduce 
small parameter p expressing the relation between characteristic longitudinal 
(denoted by A) and transverse scales, for simplicity the same for both transversal 
directions x and z. The equivalent system in the dimensionless variables : 

J?*,lc*,p*,p*,t*:J? = c 7 v^p' = c 2 pop^p > = pop^ 

_^ (-3) 
x = (Xxt/yfjl, \y* : Xz^/yfji), t = A£*/c, 



(unperturbed values marked by index zero, perturbed ones are primed, c = 

\Jl^ is adiabatic sound velocity), looks as follows (asterisks for dimensionless 

variables are omitted here and everywhere later): 

d 



dt 



1p + Dip = (f + (fitv: 



where ip is a column of the dimensionless perturbations 

T 

V X Vy V z P 



1> 



L is a linear matrix operator : 
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■8\n 



o 2 



^1 y/V- 1 dxdy 

-b l u- d2 



dxdz 

d 2 



dydz 



yjjid/dx 
d/dy 



~^i\/^d d zd y 
d/dy 



'dxdz 

JpZd/dx 

y/Jld/dx d/dy 
with dimensionless parameters originated by thermal conductivity and viscosity 



V 



AVfe - S 2 A Jpd/dz 
Jpd/dz —6% A 
^fjld/dz 



\ 



/ 



(•4) 



(.5) 



.6) 



(C + n/s) - 

5 X = 



51 = 



X 



Pqc\ 1 1 pocX' 2 pocXCy 1 z pqcXC v ^/ 
There introduced also dimensionless operators V, A : V = 

y/jid/dx d/dy y/jid/dz), A = pd 2 /dx 2 + d 2 /dy 2 + pd 2 /dz 2 . The right-hand 
side of Eq. Ofl consists of two quadratic columns, the first <p which does not 



X 



depend on thermoviscous effects 

/ , 



V 



vV)v x + y/ppdp/dx 
-{vS7)v y + pdp/dy 
vV)v z + sfppdp/dz 
-jp(Vv) — (vV)p 
-p(W) - (vV)p 



\ 



(•7) 



/ 



and the second one ip tv appearing in the thermoviscous flow: 



/ 



-5fpAv x - 5\pUVv 



\ 



dx 



-5\pAv y - 8\pf y (Vv) 



-5\pAv z - 6lp£ z (Vv) 

( (7 - 1) ((si - m (w) 2 + f + n - f^ny 



Linear flow is defined by the linearized version of the system of EqsOl) 

J^ + 2^ = 0. (.9) 

For linear flows, a solution may be found as a sum of planar waves: v x = 
v x { k ) exp(iut — i k 1?), .... with wave vector k = [k Xl k y , k z ). In the space of 
Fourier transforms (marked by tilde), — ik x means d/dx, (—iky)^ 1 represents 
f dy, iu means d/dt, etc. System (C3 ) yields in the five roots of dispersion 
relation: 



ui = Q + u 2 = -£l + i(3ti z /2, cu 3 = -z^Q 2 , u 4 = u 5 = iS'ftt 2 , (.10) 



where j3 = 5\ + 5\ + 5\ + 5%, Q = k y + 



These five frequencies relate 



to three branches of possible motions in fluid, two acoustic modes, the entropy 
mode and two vortical modes. For the real substances, (3 > 0,£ 2 > and 
8\ < 0, that provides correct signs of imaginary parts of all frequencies. Modes 
as eigenvectors of a linear problem in the space of Fourier transforms look: 
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v x i ( k x i ky , k z ) 
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.11) 



As basic variable (which all other perturbations are expressed through) the 
Fourier transform of density perturbation is chosen for the first three modes. 
Vorticites 04, 05 are chosen for the last two because the vortical modes possess no 
density and pressure perturbations. All calculations of the modes and projectors 
have accuracy up to the terms of order fi, (3. Specific features of the modes 
follow from Eqs d.lll ): both vorticity modes keep unperturbed density, entropy 
mode is isobaric, and so on. Any field of the linear flow is sought as a sum of 
independent modes. To calculate projectors that decompose a concrete mode 
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from the overall field ip, matrix M has to be defined: 



M = ( Pi P2 P3 <Jh) T = V> = 



(.12) 



as well as the inverse matrix M 1 
perturbations into specific modes 



Projectors that decompose the overall 



Piifj = if)i,i = 1, ..5 



(.13) 



may be calculated as a product of a column with number i of M and row 
number i of the inverse matrix M^ 1 : 



P i = M i -M i ,z = l, ..5. 



(.14) 



Calculations provided with accuracy of order fj.,/3 yields in projectors as 
follows: 



Pl,2 = 



\F&2k v 







2 ) 



A = 



%b\k y 



^0 






^0 i(8\ + 8l)ky -1 







\ 





1 



/ 



.15) 



9 



Pa 












II 



h 2 





V^feOO 

kl 



\ 



\ 



i — k 


k 2 













































Pk 



7 















V 



7 




Matrix projectors satisfy common properties of orthogonal projectors: 



yP i = I,P i -P n = if i^n,Pi'Pi = Pi if i = n,{i,n=l,..5), (.16) 



i=l 



where /, are unit and zero matrices. The inverse transformation of formulae 
( 1.151) to the (~x,t) space may be easily undertaken. 

Since modes in linear flow are exactly decomposed by projectors, the dynamic 
equation for every mode may be originated from the system Oj ) by acting of 
the corresponding projector: P (J^ + Lift) = 0. The famous dynamic linear 
equations for acoustic pressure is as follows: 



9pi,2 ± dpi,2 j_ At 



= 0. 



:-17) 



dt dy 2 J * 2 dy 2 

Equations for vortical and entropy modes are ordinary equations of thermal 
conductivity and are well-known also and may be found in many sources [16 ] . 
Acting by projectors at the original system of Eqs Ol) with non-zero nonlinear 
part (that essentially depends on all modes) results in coupled equations for 
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interacting modes: P (J^ip + Lif) = P {ip\ + (pu v ) ■ When one of the modes 
is dominant and role of all other ones in quadratic source is ignored, acting 
of the corresponding projector yields in nonlinear dynamic equations, the fa- 
mous Earnshow one for the plane flow (// — ► 0), the Khokhlov-Zabolotskaya- 
Kuznetsov (KZK) equation for acoustic pressure in quasi-plane viscous flow as 
follows: 

and other equations. 



STREAMING CASED BY ANY ACOUSTIC SOURCE. RADIATION FORCE OF 
STREAMING. 

The streaming as secondary flow induced by losses of acoustic momentum 
and its partial nonlinear transfer to the momentum of vortical mode is best 
observable when nonlinear effects may store with time. Most experiments deal 
with quasi-periodic waves or wave packets of many periods. That gives rise 
to the theoretical basis concerning namely to this kind of source. An area of 
projecting is much more extended that streaming since all possible interactions 
of modes may be calculated. From the point of view of nonlinear interactions, 
streaming is a particular case of dominant acoustic source inducing the vortical 
flow. The governing evolution equation follows when one acts by sum of projec- 
tors P4 + P5 at the system of Eqs(ClI). Only acoustic terms related to the first 
mode (rightwards progressive) should be left at the both quadratic columns . 
Physically, that corresponds to the initial stage of nonlinear interaction when 
the acoustic mode is dominant and the inverse influence of the streaming may 
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be ignored. Also, the generation of heat is ignored. A dynamic equation for 

transversal component of velocity V x looks 

r)V r) 2 V 

-W-^W = - (v ^ )Vt + Fix ' m 

where F\ x is a transverse quadratic acoustic source (radiation force) caused by 
the rightwards progressive acoustic mode. Dropping calculations, the result is : 

ix \^q x Qy 5^1 Q x Qy q x J dy{ Qy)^j (20) 

The radiation force appears in the viscous flow only. As expected, all quadratic 
nonlinear terms in the acoustic source of order (3 {) give zero in sum and therefore 
attenuation is necessary for acoustic streaming to exist. Acoustic pressure p\ 
should be a solution of KZK equation. Projectors P4P5 and their sum yield in 
the vortical radiation force automatically due to their origin. 

Let us discuss the result to demonstrate the difference between a periodic 
and pulsed acoustic sources. 



1. In a case of a pulsed ultrasound as in the papers [LL9J], [20] with cl Cell" - 
rier frequency coo the solution of the equation ( 1.191 ) with the right-hand side 
of Eq. ( 1.201 ) contains relatively small (non-growing) part, oscillating with the 
frequency 2loq. It is known, that there is also a term growing proportional to 
time, that arises from constant average of a squared oscillating functions that 
is evident from the identity Sin 2 (tJot) = (1 — Cos (2cuot))/2. If one take a wave- 
train of a rectangular form, it gives the input proportional to Pqt/2; in the case 
of a triangular pulse one gets P 2 r 2 /4 (both per pulse). The second result is 
obviously larger if r > 27t/cjq. This simple example explains why the streaming 
amplitude depends on pulse structure of the source. Experimental data prov- 
ing the larger velocity of streaming caused by pulsed mode are discussed by 
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Starritt, Duck, and Humphrey 

2. It is very important that now the averaged force (Fi x ) contains the addi- 
tional terms that cancel in the case of a periodic sound field. The integration 
by parts highlights the point: 

-y+Ay 



ry+^y 

/ ( m y+p 2 y )dy = o.5(p 2 ) y \y +A y, (.21; 



the right-hand side is obviously zero for a spatially periodic function with pe- 
riod Ay standing for the acoustic pressure. This shows non-equivalence of the 
expressions for the force obtained by averaging procedure by different authors 
delivered in different stages of the streaming equation derivation. In other 
words if one adds such expression to the force, the mean value of it survives in 

^periodicity the term contributes. 



a periodic wave case while for any deviation o 

Radiation force given by Gusev, Rudenko |l£f in the special case of strongly 
attenuated quasi-periodic acoustic beam with weak diffraction is given by 

= v^P^(0 2 /2) e W (-Py) (-22) 

rewritten in our variables. It is value averaged over integer number of sound 
periods and relates to acoustic pressure in the rightwards beam as follows: 

P 

pi(x, y, t) = P 6(x) exp(--y) sin(t - y). (.23) 

Calculating based on formula ( 1. 201) results in the following: 

F lx - $ lx = 5^ sin(2(t - y)) + 0(p 2 ). (.24) 

Averaging over integral number of period of the sound wave gives (F\ x ) = $1^ 
in the leading order. 
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NUMERICAL EXAMPLES AND DISCUSSIONS. 



In the frames of report by Gusev, Rudenko [18], a strongly attenuated non- 
diffracting acoustic beam in unbounded volume is considered. The calculations 
give values of frequencies above 1MHz in air in the case of radiator with a ra- 
dius of order 0.1m. That considerably simplifies the further calculations since 
a source should be thought as a solution of the KZK equation and till now does 
not have an appropriate mathematical basis. Even the Kholhlov-Zabotskaya 
(KZ) equation suitable for flows without attenuation is extremely difficult for 
analytical solution: evaluating of flow in the paraxial region needs many inter- 
mediate calculations | . In spite of general formula on radiation force given 
by Eq. fl.20p . any simple solution is desired for illustrations though the stream- 
ing caused by any other wave satisfying the KZK equation may be treated by 
this formula. As an appropriate acoustic source satisfying the limitations of 



paper by Gusev, Rudenko 
as follows |2|: 



let us take a mono-polar two-dimensional wave 



I ft? I 2^ e *P (~^) ,,,, 

where £ = (3y, r = t — y, e = ^±1. A self-similar solution like ( 1.251) possesses 
correlated values of characteristic scale ad amplitude in contrast to periodic 
perturbations where effective Mach number and wave length may be sought 
independent. Constant C is responsible for the shape of the curve: large C 
provides a curve close to the Gauss one. To calculate a transverse velocity 
Vix, the first, corresponding force F\ x should be calculated, and the second, the 
momentum equation ( 1.191) has to be solved. 

The formula of this chapter apply to the three-dimensional flow and are also 

14 



suitable in the two-dimensional or plane one. Generally, two parameters /i x ^z 
instead of one \x responsible for the geometry of flow may be introduced from 
the very beginning. To consider two-dimensional geometry, the flow over z-axis 
supposed to be uniform. To calculate the field of velocity accordingly to this 
last equation, let us do some remarks. It is important that the hydrodynamic 
nonlinearity which is necessary in periodic sources (otherwise the streaming 
grows infinitely [ijj j3| 2|) may be not considered in the pulse dynamics, 
because of absence in this case the storage of nonlinear effects. Amplitudes of 
the generated streaming is so small that the quadratic term such as this one 
may be omitted. Without the mentioned term, Eq. ( 1.191) is a simple equation of 
the heat conductivity with viscous coefficient cased by shear viscosity 5\ with 
an acoustic source in the right-hand side. 

In order to simplify calculations, the viscous term in Eq.( l.l9l) will be ignored. 
The same limit relating to low shear viscosity was considered in the referred 



paper by Gusev, Rudenko 



181 ] . Note that viscous terms grow significantly near 



the discontinuity of the profile of V x , that never happens to a smooth profile 
caused by pressure source given by Eq.( l.25l) . To calculate the temporal field of 
velocity at a set of co-ordinates y, the integration of the radiation force F\ x has 
been occurred: 

V x = [ F lx dt. (.26) 
J o 

For any co-ordinate y, ignoring of viscous term results in stationary velocity 
after the source passed. The reason is the constant result of integration of the 
source from minus infinity till plus infinity. Indeed, there is an attenuation 
in the very Eq. ( 1.191 '). not only in the attenuated acoustic source, so for the 
more detail calculations it should be involved. The radiation force also needs 



15 



integration over y, the constant of integration is chosen to be zero when y — > oo. 
Transversal derivatives in formula for the radiation force and V x are simple when 
diffraction is ignored: 

d J^MA = -2 X eM-x\ (-27) 

so all calculations of velocity in this paper refer to the point x = ^ where the 
radiation force achieves maximum. Constants of Eq. ( 1.251) used in calculations 
are: C = 2, e = 1.2, (3 = 0.1 

Figure 1 shows the temporal development of velocity V x of streaming divided 
by J~jif3 (one may get dimension values by multiplying by c) at different points 
y:l, 2, 3,.. 10; x = y/2/2. The stationary field of after passing the source is 
presented by Fig. 2. It is evident that velocity of streaming tends to a constant 
level with time passing in contrast to the streaming caused by periodic source 
when hydrodynamic nonlinearity ignored. 

The next, the streamlines may be plotted due to the calculated V x and the 
known relation for the vector field: V V = 0. The following integration with the 
proper constant should be occurred over y to get a field of longitudinal velocity 
V y . Finally, the calculated streamlines at t = 3,t = 5 are shown at the figures 
3, 4(a) both with the dimensionless pressure of acoustic source correspondingly 
(Fig. 3, 4(b)). The symmetrical lines of the curves are x = 0.5 at the upper 
half-space and x = —0.5 at the lower one due to the chosen shape of the source 
given by exp(—x 2 ). 

Some calculations of dimensional parameters of flow seems to be useful under 



limitations of report by Gusev, Rudenko |l8j ]. Quantitative evaluation relates 
to air. For example, a characteristic length of source is A ~ 5 • 10 _5 m, that gives 
a dimensionless (3 = 4 • 10~ 6 m/A ~ 0.1. A radius of transducer is R = 0.1m. 
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Pressure of acoustic source and therefore streaming depends on small param- 
eter P; the second one , \x appears as a multiplier in the expression for the 
radiation force. If hydrodynamic nonlinearity and shear viscosity ignored, V x 
is simply proportional to yfjl = A/R. Accordingly to the accepted values of A 
and R, y/Ji = 5 • 10~ 4 . Calculations show that for effective Mach numbers of 
source about 4- 10~ 2 velocity of streaming achieves stationary level 5- 10~ 4 m/sec 
(see figures 1,2, y=3). Meaning a single pulse as acoustic source large veloc- 
ities of streaming are hardly expected while there is no storage of nonlinear 
transport of acoustic momentum over many periods. Nevertheless it looks not 
extremely small in comparison to the measured values of streaming that vary 
from 10~ 3 m/sec till lm/sec [8]. The sensitivity of modern technique enabled 



streaming velocities down to 10 m/ 'sec 20] 



CONCLUSIONS 

The most important result of the projecting is to get space and temporal 
structure of any mode independently of the type of source. The basic idea 
of the projecting starts from separating modes accordingly to their specific 
properties in the weakly nonlinear flow. These basic motions (or eigenvectors of 
the linear flow) should be defined at first. It may be proceeded in the algorithmic 
way for both homogeneous and inhomogeneous backgrounds ( see paper on 
interacting modes in bubbly liquid |nj|), media affected by external forces and 
the real geometry of the flow. The definition of modes is unique, determined 
by the linearized differential conservation equations only. The flow in bounded 
volumes or relating to problems that need special conditions has to be sought as 
a superposition of specific modes corresponding to the concrete problem. Any 
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mode is distinguished by a correspondent projector at any time. Projectors 
apply to arbitrary source and type of initial perturbation. 

The next, the interaction of modes in weakly nonlinear flow yields in coupled 
nonlinear equations for the modes that may be solved approximately |H|, |l4|. 
Moreover, the evolution equations may be corrected up to the higher order 
nonlinear terms due to increasing influence of the other generated modes. In 
the present paper, rightwards acoustic mode is sought as a dominant one which 
gives rise to the vortical one. Therefore, the governing equation ( 1.191 ) involves 
a pure quadratic acoustic source. For more advanced flow, the growing role 
of the non-dominative modes due to nonlinear interactions may be accounted 
as well. The similar calculations were undertaken while one-dimensional flow 



studied [111]. 

Since the method applies to flow with any initial conditions (including non- 
acoustic) and does not need quasi-periodic sources, it is useful for investigation 
of some special problems like streaming caused by the non-periodic mono-polar 
acoustic source. The possibilities of analytical methods in the study of stream- 
ing are superior over experimental and purely numerical investigations in the 
view of complexity of the whole phenomenon. In this paper, the general for- 
mulae on interacting modes are presented, and the particular case of streaming 
caused by mono-polar source is discussed and illustrated by numerical calcula- 
tions. A special meaning of the quasi-periodic and particulary pulse ultrasound 
is its importance in medicine, in delicate exploring of the parameters of fluid 
as well as importance of the secondary modes following the source, e.g. in the 
sonochemistry (21] . 
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Figure Captures 

FIG. 1. Dimensionless transversal velocity via time at different longitu- 
dinal points y:l,2,3,..10 (from the lowest to the upper curve), x = y/2/2. 

FIG. 2. Stationary dimensionless transversal velocity via longitudinal 
co-ordinate y at x = y/2/2. 

FIG. 3. a)Streamlines in the plane (x,y), t=3; x-transversal and y - 
longitudinal co-ordinates, b) dimensionless pressure of acoustic source via y at 
t=3, x = y/2/2. 

FIG. 4. a)Streamlines in the plane (x,y), t=5; x-transversal and y - 
longitudinal co-ordinates, b) dimensionless pressure of acoustic source via y at 
t=5, x = y/2/2. 
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